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Abstract. In this paper we study the numerical passage from the spatially homogeneous 
Boltzmann equation without cut-off to the Fokker-Planck-Landau equation in the so- 
called grazing collision limit. To this aim we derive a Fourier spectral method for the non 
cut-off Boltzmann equation in the spirit of | 21 | , | 23 | . We show that the kernel modes that 
define the spectral method have the correct grazing collision limit providing a consistent 
spectral method for the limiting Fokker-Planck-Landau equation. In particular, for small 
values of the scattering angle, we derive an approximate formula for the kernel modes of 
the non cut-off Boltzmann equation which, similarly to the Fokker-Planck-Landau case, 
can be computed with a fast algorithm. The uniform spectral accuracy of the method 
with respect to the grazing collision parameter is also proved. 

Key words Spectral methods, Boltzmann equation, cut-off assumption, Fokker-Planck- 
Landau equation, grazing collision limit. 



1. Introduction 

In this paper we are concerned with numerical resolution methods for collisional equations 
arising both from kinetic theory of rarefied gases and plasma physics. These equations take 
the form 

(1) d t f + vV x f = Q(f,f), t>0, x,veR 3 

where f(x,v,t) is assumed to be a nonnegative function that represents the density of 
particles in position x at time t with velocity v. In ([T]), Q(f,f) is a quadratic integral 
operator describing the collisions of particles, whose form we shall make explicit below. 
Here we will mainly focus our attention to the approximation of Q(f, /). 

As it is well known in the physics literature, the Fokker-Planck-Landau (FPL) equation 
is used to describe the binary collisions, occurring in a plasma, between charged particles 
and can be written as 

(2) d t f + v-V x f = Q L (fJ), 

(3) Ql(IJ) = Vv( [ dv*a(v-v*)[f*V v f-f(V v f), 
where a is the matrix-valued function of the form 

aij(z)=n*\)(k3-T&), 
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and we use the standard notations 0* = 4>(x, i>*, t). The kernel "J depends on the interaction 
between particles. A typical choice is, up to a multiplicative constant 

*(|z|) = |zp+ 2 , -3< 7 <1. 

This operator was introduced as an approximation of the Boltzmann collision operator 
in the case of Coulomb interaction (i.e. 7 = —3) [TB]. The motivation is that the Boltz- 
mann operator is meaningless in the case of a Coulomb interaction. In fact, for a Coulomb 
interaction, the effect of those collisions that are grazing, i.e. collisions which result in an 
infinitesimal angle deflection of the particles trajectories, prevails over the effect of other 
collisions. 

To clarify this idea, let us consider the Boltzmann equation in the form [5], [5] 

(4) d t f + v-V x f = Q B (fJ), 

(5) Q B (fJ)= [ dvJ d*B(v-v.,9)(f'fi-ff.), 

Jr 3 js 2 

where the deflection angle 9 <E [0,7r/2], is such that cos# = (v — u*) • a/\v — The unit 
vector cr parameterizes the set of all kinematically possible (i.e., those conserving energy 
and momentum) post-collisional velocities v' and v'^ by 

, v + u* \v — vJ 
= + —2—* 

\v - v»| 

-a. 



2 2 

The relative probability of these outcoming velocities depends on the nature of the interac- 
tion between the molecules, and this is taken into account in the kernel B. For molecules 
which interact through an l/r s force law, where r is the distance between interacting par- 
ticles, physical arguments show that the natural choice for the kernel B is 

v — v 

B(v — v*,0) = \v — v*| 7 &(cos0), cos# = (- -t,c), 

\v-v*\ 

where 7 = (s — 5)/(s — 1), and sin 6b(cos6) is a smooth function except near 6 — 0, where it 
presents a (nonintegrable) singularity of order (s + l)/(s — 1). This singularity corresponds 
to the grazing collisions and when these collisions prevail (in a way we shall make precise in 
the next section) solutions to the Boltzmann equation converge towards solution of the FPL 
equation. In recent years, a noticeable amount of theoretical results [2] , [TD] , [T2] , |30 j .|31 j 
concerned with the limiting process involved, called the grazing collisions asymptotics, has 
been obtained. 

On the other hand, the usual numerical resolution methods for the Boltzmann equa- 
tion can not be used in this context, exactly in reason of the singularity in the kernel B. 
Classically, the standard strategy for numerical applications is to cut-off the small angle 
collisions so that B becomes integrable. In fact, under this simplification, it is possible to 
split Qb{J, /) into a gain term Qg(f, /) and a loss term Qg{f, f, ) according to 

QB(/,/) = Q^(/,/)-Qi(/,/) ! 

which makes the problem treatable both by direct simulation Monte Carlo methods [5].|17].[2"0] 
or fully deterministic methods [1].|26|. In particular, as the cut-off parameter becomes 
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smaller, the computational cost of these methods increases dramatically so that they are 
unable to treat in practice the non cut-off problem. Recently a study of the behavior of 
Monte Carlo methods for a simplified Boltzmann model close to the non cut-off case has 
been presented in [13J . 

To the best of our knowledge, the method we have studied in the present paper is the 
first method which is able to treat numerically the Boltzmann collision operator in the non 
cut-off case. 

In very recent years, there has been a new approach, based on a Fourier spectral method, 
to the numerical approximation of the Boltzmann and the FPL operators [21] , [23] , [24] . The 
advantage of this method is that the whole structure of the collisional operator (in practice 
all the informations characterizing the kinetic equation) is contained into a series of kernel 
modes. In addition, in the FPL case, thanks to the structure of these kernel modes, the 
resulting algorithm can be evaluated with a strong reduction in the computational cost. 
Other numerical methods, that avoid the quadratic complexity of the FPL operator, have 
been introduced recently [5].|ll j .[19 ] . 

Here, we consider this spectral method in the case of the Boltzmann collisional operator 
([5]), without the cut-off assumption, by showing that the corresponding kernel modes are 
well-defined for any reasonable kernel with singularity. The main result here is that these 
kernel modes converge (in the grazing limit) to the corresponding kernel modes of the FPL 
equation and that the spectral method is uniformly accurate with respect to the grazing 
collision parameter. In addition, we derive an approximate representation for the kernel 
modes of the non cut-off Boltzmann equation which is valid for small values of the deflection 
angle and that, similarly to the FPL case, can be computed through a fast algorithm. 
Numerical results based on this representation are actually under development and will be 
presented elsewhere. 

Summarizing the method has the following properties: 

• Well-defined for the non cut-off Boltzmann equation. 

• Correct numerical grazing collision limit. 

• Uniform spectral accuracy with respect to the grazing collision parameter. 

• Strong reduction of the computational cost close to the grazing limit. 

For the sake of completeness, we mention here that the idea of using the numerical 
asymptotic limit as a guide for the development of numerical methods is related to some 
recent works in the context of the fluid- limit of the Boltzmann equation (see [6], [15] and the 
references therein). 

We briefly discuss the organization of the paper. Section 2 is devoted to recall some 
mathematical and physical properties of the equations and to present the precise formulation 
of the grazing collision limit. Section 3 deals with the derivation of the spectral projection 
of the non cut-off Boltzmann equation. The grazing collision limit of the kernel modes is 
contained in Section 4. Finally, Section 5 is devoted to a proof of the uniform spectral 
accuracy of the method. Some final considerations are contained in Section 6. 

2. MATHEMATICAL AND PHYSICAL BACKGROUND 

In this section, we briefly recall the mathematical and physical background both of the 
FPL and Boltzmann equations. In the rest of the paper we will assume a kernel B with a 
non integrable singularity, as described in the introduction, and we will consider the space 
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homogeneous equations. It is well-known, in fact, that by the standard splitting algorithm 
we may restrict to the space homogeneous case. 

From the physical point of view, the most important quantity associated to the angular 
cross-section 6(cos 9) is the so-called (angular) cross-section for momentum transfer, which 
is (up to a multiplicative constant) 

r r*/2 
A= / b(k ■ cr)(l - k ■ a)da = 2tt / &(cos 9) (1 - cos 9) sin 6d6, 
Js N - x Jo 

or more generally the function 

(6) A{q) = 2Tr B{q,0){l - cos 0) sin 9 dB. 

Jo 

It is shown in [30 that for — 3 < 7 < 1, the finiteness of A is exactly what is needed to 
develop a mathematical existence theory. 

If s = 2 (Coulomb interaction), then A is a divergent integral as 9 — > 0, which prevents 
to define Qs(/, /) in any reasonable sense (see [29j Part I, Appendix A]). Using heuristic 
arguments, Landau showed that the principal part of ([5|) under the truncation (9 > e > 0) 
is proportional to Ql(J, /) (see [55] for a rigorous variant of Landau's argument). 

In view of these remarks, the following setting will appear to be natural from both the 
mathematical and physical points of view. We use the standard notations 

L\ = U £ L\R 3 ), ^ |/(w)|(l + \v\ s )dv < ooj , 

L\ogL= j/GL 1 ^ 3 ), jfj/(«)|bg(l + \f(v)\) dv< ooj. 

Moreover, if / is a nonnegative function, we shall call J f dv its mass, J fv dv its momentum, 
/ f{\ v \ 2 /2)dv its energy, and J flogfdv its entropy. 

Definition 1. We shall say that the family of kernels {B e ) e>0 given by 

(7) B e (q,6) = \q\n e (co S e) 
(—3 < 7 < 1) concentrates on grazing collisions if 

a) A £ = 2ir 6 e (cos 0)(1 - cos 9) sin 9 d9 — > A < 00, 

^0 

b) V6»i > 0, b e (cos9) — ^0 uniformly on (9 > 9^. 

In fact, we could allow much more general assumptions on B e . For instance, 
(i) = 

J B E (q, 9) sin 3 9 d9<C^ + \q\ 



(ii): 

P Tr/2 



2n\q\ 2 / B e (q,6)(l -cos 9) sin 9 d9 A (|g| 
Jo 



in< o f (M 3 \W), 
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(iii): 

, ,7T/2 

Vi?>0, / |g| 2 / B e (q,0) sin 4 0d8dq — ► 0. 

J\q\<R JO 

And then, the function -Ao ( I Q I ) would replace the function Ao|<7| 2+7 in all that follows. 
The main results of [3D] can then be stated as follows, 

Theorem 1. Let /q be an initial nonnegative datum with finite mass, energy and entropy. If 
7 > 0, assume in addition that fo £ L^ +s for some S > 0. Then, for all e > there exists a 
solution f £ of the Boltzmann equation with kernel B e , such that / e (0, •) = fo- Moreover, as 
e — > 0, up to extraction of a subsequence, for all T > f e converges weakly in L 1 ([0,T];K 3 ) 
to a solution f of the FPL equation with kernel ^(1^1) = ^A|z| 7+2 . 

Remark 1. The exact value of the cross-section given in the previous theorem, comes out 
in a clear way from the analysis of the numerical grazing collision limit we will develop in 
Section 4. 

We do not make precise here the meaning of weak solutions, which is not the usual one 
for 7 < —2 (see [30] for details). 
From now on, we use the notation 

C(6>) = b(cos6) sinO, 

in which the Jacobian arising from the spherical variables is taken into account, and we 
recall that this function has typically a sing ularity at 9 = like 0- ( - 1+u \ v = 2/(s + 1) for 
inverse s-forces. Typical families (£ e ) are given by 

U0) = T^—A e ) 1 ^ f°r S = 2, 
lne^ 1 

< 

a0)sin 2 ^ = lsin 2 ^C(^ for S >2. 

We emphasize that from the physical point of view, the case 7 = —3 is the most interesting 
(see also |27j . in which a somewhat different viewpoint on the FPL equation is given). 

Now, let us give the state of the art concerning the mathematical theory of the spatially 
homogeneous FPL equation. The first works on the subject were done by Arsen'ev and 
Buryak [2 . Their results were considerably improved and extended in the recent work [14] , 
from which we extract the following theorem, which gives a quite satisfactory picture of the 
case 7 > 



Theorem 2. Let 7 > 0, and let fo be a nonnegative initial datum, /o € f or some 

6 > 0. Then, there exists a solution f(t,v) of the FPL equation, continuous in time, with 
/(0,-) = fo, satisfying the conservation of mass, energy, and the decrease of entropy, and 
which is very smooth for all positive time. More precisely, for all t > to > 0, and k, s > 0, 
one has f £ H k C\L\, where H k denotes the standard Sobolev space of order k, and the norm 
of f in these spaces depends only on fo and to- In addition, if J |/o(w)| 2 (l + \v\ s ) < 00, with 
s > 15+57, then f is unique (in the class of weak solutions which may let the energy decrease) 
and for all time t > 0, there exist positive constants A t and K t such that f(t, v) > A t e~ Kt ^ . 
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In the case 7 = 0, it is sufficient to require that / G L\ to get both existence and 
uniqueness of solutions [31] . 

As regards other qualitative features of solutions, a detailed study of the asymptotic 
behavior as t — > 00 was done in [14], which led to the following results 

Theorem 3. Assume 7 > 0, and let f be a (smooth) solution of the FPL equation with 
kernel |z| 7+2 . Without loss of generality, assume that f has mass 1, momentum and 
energy 3/2. Then, 

||/(t,-)-M|| i x<C7(/o)*- 1 /' 1f , 
where C(/o) is explicit and depends only on the entropy of fo, and 

Miv) = ww- 

Moreover, f — M goes to in all (weighted) Sobolev norms. 

Again, in the case 7 = the results are simpler, in that the trend to equilibrium is 
exponentially fast, with a bound of the form Cee"^ 2-5 -**, for all e > 0. 

As a consequence of these results, it was shown in 14 that if 7 > 0, and if the initial datum 
satisfies the conditions for uniqueness of the solution, then this solution is stable, globally 
in time, say in L 1 norm, with respect to small perturbations of both the initial datum (in a 
suitable weighted L 2 space) and the kernel, in the sense that one allows perturbations of the 
form ^(|^|) = |z| 7+2 (l + 7y(|z|)), with r\ small in C 2 norm. From the practical viewpoint, 
this means that small errors in computations arc allowed. 

We conclude this presentation by noting that in the case 7 < 0, the mathematical theory 
is still at an early stage (only existence of global weak solutions is known, and the study of 
trend to equilibrium is still in progress). 

3. Spectral projection of the non cut-off Boltzmann equation 

In this section we derive the Fourier spectral method for the non cut-off Boltzmann 
equation following [2"T ] . [2"3" ] . 

A simple change of variables permits to write the Boltzmann collision operator Qs(/, /) 
in the form 

(8) Qfl(/,/) = / dq [ daB( q ,e)[f(v + q+)f(v + q-)-f(v)f(v + q)}, 

JR3 JS 2 

where q — u* — v and the vectors q + and q~ that parameterize the post-collisional velocities 
are given by 

(9) Q + = \(q+\qW), <T = \{q - WW). 

We point out that the possibility to integrate the collision operator over the relative velocity 
is essential in the derivation of the method. Other kinetic equations with a similar structure 
may be treated in the same way. 

First we need to reduce the problem to a bounded domain in d. To this aim we observe 
that if a distribution function / has compact support, Supp(/(t>)) C B(0, R), where B(0, R) 
is the ball of radius R centered in the origin, then by conservation of energy and momentum 
the collision operator satisfies the following [21] 
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Proposition 4. Let Supp(/(t>)) C 23(0, R) then 

QB(fJ)(v)= [ dq [ daB(\q\,6)[f(v + q + )f(v + q-)-f(v)f(v + q)], 
JB(0,2R) Js 2 

with v + q+,v + q-,v + qe 23(0, (2 + \/2)R). 

Thus to develop a spectral method we will consider the distribution function f(v) re- 
stricted on the cube [-T, T] 3 with T > (2 + V2)R, assuming f(v) = on [-T, T] 3 \ 23(0, i?), 
and extend it by periodicity to a periodic function on [— T, T] 3 . 

As observed in [23] using the periodicity of the function it is enough to take T > (3 + 
\/2)R/2 to prevent intersections of the regions where / is different from zero. 

To simplify the notation let us take T = tt and hence R — \ir with A = 2/(3 + V2). 
Hereafter, we use just one index to denote the three-dimensional sums with respect to the 
vector k — (ki,k2,k 3 ) £ 2Z 3 . 

The approximate function is represented as the truncated Fourier series 

JV 

(10) M«)= ]T he ik - v , 

k=-N 




The Fourier-Galerkin method [7].|16] is obtained by considering the projection of the homo- 
geneous Boltzmann equation on the space of trigonometric polynomials of degree < N. 
Hence, we have 

(12) / [d t f N -Q B (fNjN)]e- ik - v dv = 0, k = -N,...,N. 

J[-n,ir]3 

By substituting expression (jTu| in (IT21 we get a set of ordinary differential equations for the 
Fourier coefficients 

N 

(13) d t f k = Y, flfmP(l,m), k = -N,...,N, 

l,m=-N 

l + m = k 

where the Boltzmann kernel modes 0(1, m) are given by 

(14) /3(/,m)= / dq I daB(q,d) 

JB(0,2\tt) Js 2 

Note that, as a consequence of the properties of trigonometric polynomials, di4]) is a 
scalar quantity completely independent on the argument v, depending just on the particular 
kernel structure. This property is strictly related to the use of a Fourier spectral method. 
Other spectral methods may be developed, however they do not lead to this fundamental 
simplification. 

In practice all the informations characterizing the kinetic equation are now contained in 
the kernel modes. Clearly, these quantities can be computed in advance and then stored in 
a multidimensional matrix. 



il-q -\-im-q 



a iq-m 
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Finally, the integral in (|14[) can be written in the form 

$(l,m)= / / B(p,a-ui)e ip ^ l+m ^ e -w(.\l-™W-i>) _ e -i P (\i-m\w) pdpdoJ da, 

Jo Js 2 Js 2 L J 

where q — puj, (I + to) — \l + m\r\ and (/ — to) = |/ — m\p. 

It is possible to show that the previous expression is invariant under any unitary transform 
R (rotation or symmetry) of the unit vectors 77 and \i [21] . As a consequence the kernel modes 
are functions only of \l — m\, \l + m\ and of the angle between rj and p. In particular we have 

(15) p(l,m) = F(\l + m\,\l-rn\,r)-iJ,) = F(\l-m\,\l + m\,r}-ti). 



Since from (TT4"]) it follows easily [3(1, to) = p(— I, — m), from ([15]) we get also 0(1, m) — 
j3(—l, —to) and hence the kernel modes are real functions. 

Obviously, these properties are useful to reduce the storage requirements of the method. 
We refer to [53] for a more detailed discussion on this topic. 

Remark 2. 

i) The evaluation of (fT"3"]) requires exactly 0(N 6 ) operations. We emphasize that the 
usual cost for a method based on ^V 3 parameters for / in the velocity space is 
0(n a N 6 ) where n a is the numbers of angles discretizations. Thus the straightfor- 
ward evaluation of ()13[) it is less expensive than a usual discrete-velocity algorithm. 

ii) In the VHS case, B(q,9) = C a \q\ a , the dependence on the scattering angle disap- 
pears and ([M)) reduces to a one-dimensional integral [23] 

/> Xtt 

(3(1, to) = C x , a / r 2+a [SincflJ + m\r) Sinc(|Z - m\r) - Sinc(2|m|r)] dr, 
Jo 

where C\ yCC — (47r) 2 (8) a C Q . In addition for integer values of a < —3 (like in the 
case of Maxwell molecules or hard spheres) the previous integral can be computed 
explicitly [23] given rise to very simple and easy implementable approximation for- 
mulas. 

4. The grazing collision limit of the kernel modes 

From the previous section it is clear that in the grazing collision limit it would be desirable 
that the previous spectral projection gives a consistent spectral method for the limit FPL 
equation. 

To this aim first we compute the grazing collision limit of the Boltzmann kernel modes 
and then we verify directly on the FPL equation that the asymptotic limit of the Boltzmann 
kernel modes coincides with the kernel modes of the FPL equation. 

4.1. Computation of the grazing collision limit. We shall compute the limit of the 
Boltzmann kernel modes, as the parameter e defining the asymptotics of grazing collision 
goes to 0. 

Thus we start from (|14p written in the form 

e * ( '~ 2 m) -(g-kk) _ 1 



(16) f3 E (l,m)= dqe lqm da B £ (q,0) 

JB(0,2Xit) Js 2 

where B e (q,9) is given by ([7]). 
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Let us introduce a system of spherical coordinates around the vector q. Hence (q/\q\,j, h) 
will be a fixed basis, and the spherical angle <j> will be defined by 

\q\a = qcosd + j\q\ sin0sin(/> + h\q\ sin0cos<^>. 

We get 



Pe(l,m)= / dqe iqm d8sin8B £ {q,8) d<j)[e l 

Jb(0,2\tt) Jo Jo 



with 

n _ m \ 

(17) Fi- m (q, 0, 4>) — i • [q(l — cos 9) + j\q\ sin 8 sin <j> + h\q\ sin 8 cos tf>]. 

To start with, we use the Taylor expansion 

(18) e Fi- m (*M) _ 1 = Ft _ m ( q: 9, 0) + l -[F^ m {q, 8, 0)] 2 + 0(\q\ 3 \l - m\ 3 8 3 ). 

By the assumptions in Definition 1 (or point (Hi) after Definition 1), the remainder can be 
neglected since as e — > we have 

/ dqe lqm \q\ 3 [ d8 sin 8 B £ (q, 8)8 3 0. 

Jb(0,2\tt) Jo 

Now we consider the contribution of the first order term in (|18p . By symmetry we have that 
i\q\ / d(f) • [j sin 8 sin0 + h sin 6* cos 4>] = 0. 



2 

Thus the first order contribution reduces to 



(19) iri [ dqe iqm {l-m)-q [ d8B £ (q, 8) sin 0(1 - cos0). 

Jb(o,2\h) Jo 

Next we consider the second order term in (|18|) . Again by simple symmetry considerations 
the terms in sin 2 0sin0cos</), sin0sin</> and sin cos disappear. Moreover also 

-~ f dqe tq - m [(l-m)-q} 2 f d8B £ (q, 8) sin 0(1 - cos 0) 2 , 

is small of order e 2 since (1 — cos$) 2 is of order 4 . 
For the remaining terms we have 

i r 27r 

--\q\ 2 sin 2 8 / d(/){[(/-m)-jfsin 2 (/.+ [(;-m)-/i] 2 cos 2 (/>} 



= ~\q\ 2 sin 2 8{[(l - m) ■ j} 2 + [(I - to) ■ h} 2 } 

= ^\q\ 2 sin 2 8[(l-m^} 2 , 

where (/ — to) = (/ — m) — ((I — m) ■ q/\q\)q/\q\. 
Thus the second order contribution is given by 

(20) -J / dqe lq - rn \q\ 2 [{l-m)^] 2 f d8B £ (q, 8) sin 3 0. 

° Jb(0,2\it) Jo 
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Finally, adding (TlT)]) and (|20l) we obtain the approximation for small e 

(21) $ e (l,m) = 2ir [ dqe^ m [ d#sin#sin 2 ^-B e (q,9) Gi- m (q,6) +0(e), 

Jb(o,2\it) Jo z 

with 

(22) Gi- m (q,e) =iq-(l-m)-^cos 2 ^\q\ 2 [(l-m) ± } 2 . 

Since by assumptions in Definition 1 (or point (ii) after Definition 1), we have as e — > 

27rsin(9(l -cos6)B e (q,6) -» Ao|gf 5 e=0 , 
in the limit we obtain the following expression for the grazing kernel modes 

(23) A>0,m) = ^/ dqe i "- m \qp{iq-(l-m)-hq\ 2 [(l-rn) ± } 2 }. 

L JB{0,2\n) 4 

Now in (|23[) we have a balance between an advection term, the one with iq ■ (I — m) and a 
diffusive term, the one with |g| 2 [(? — m)^] 2 . 

Taking into account the symmetries of the kernel modes from (|15[) we have e (l,m) — 
$ e (—l,m) that in the limit e — > gives the equivalent representation 

(24) p o{hm) = -h [ dqe^ m \qp{tq-k+hq\ 2 [k^}, 

1 JB(0,2\n) 4 



where k = I + m. 



Remark 3. From the previous results some remarks are in order. 

i) The Boltzmann kernel modes (1161) are bounded whenever the family of kernels B e 
satisfy the conditions of Definition 1. In particular, taking e = 1, we proved that 
the Boltzmann kernel modes are bounded for any kernel B with finite angular cross- 
section for momentum transfer. 

ii) The final expression (|24|) is much simpler than the starting expression (|16| . It is 
possible to show [24] that this simplification allows to compute the resulting final 
algorithm in only 0(N 3 log AT) instead of 0(N 6 ) operations (see section 4.3). 

iii) Similarly, also the approximation defined by (12HI22[) . that can be used to study 
numerically the behavior of the non cut-off Boltzmann equation when the collisions 
become grazing, can be computed with only 0(N 3 log N) operations. 

4.2. Direct derivation of the FPL kernel modes. We show here that the grazing kernel 
modes computed in the previous paragraph coincides with the FPL kernel modes that we 
can derive by direct application of the spectral method to the FPL equation. 
To start with, for any test function tp we have the identity 

Q ] lU ', f)f{v) dv = —- I dv I dv* a(v - u„) • (V - V*)//*[Vy> - V(^*] 

1 JR 3 JR 3 

where 



dv I dv* //» \ b(v - v*) ■ [Vip - Vy>»] + dij(v - v*)-[dij<p + dij(p* 



bi(v — v*) = / J dia,ij(v - v*). 
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By symmetry we get 

Ql(J, f)<p(v)dv = I dv I dv* //* [2b(v - u*) • + aij(v - v«)dij<p\ 



0«(*) = *(M) [Sii-Ttf), 



From 

we have bi(z) — — 2z ^^P (in dimension n ^ 3 the factor 2 has to be replaced by (n — 1)). 
Thus the kernel modes are obtained from 



Q L (e a - v , e irn - v )e- ik - v dv = f dqe ll - v+m - v * ^f^-{-4iq ■ k - q 2 'k : f}c i 

[-ir,7r] 3 JB(0,2Att) \Q\ 



with k x = k-k-(q/\q\)(q/\q\). 

Finally we get the expression of the FPL kernel modes 

(25) L (l,m) = -4[ dq e^^M {iq . k + h q f[ k ^ }) 

Jb{o,2\k) M 4 

which is the same as (I2U1) for 



(26) *(\ q \) = h\ q \i+*. 

Finally, we have proved the following 

Proposition 5. Let (3 e (l,m) be the kernel modes of the Boltzmann equation with kernel 
B e (q,9) — |g| 7 £ e (#), where 2irQ(9)(l — cos#) is an approximation of Aq5o = o in the sense of 
Definition 1. Then we have 

lim/3 e (/,m) = f3 L (l,m), 

e— >0 

where /3l(/,?ti) are the kernel modes of the FPL equation with kernel given by V2b}). 
Remark 4. 

i) In arbitrary dimension n^3 and with a cross-section like $(|(}i) instead of |g| 7 , the 
formula is 

*{\<l\) = J7^-M\<l\)\<l\\ 

4(n — 1) 

as shown under extreme generality in [T]. We note that the function ((9) should 
of course now be equal to 6(cos6*) sin™ -2 6, and that the constant in front of the 
definition of A, Ao should be IS*™ -2 ! instead of 2ir. 

ii) It is possible to obtain several representations of the kernel modes accordingly to 
the several possible ways in which the FPL equation can be written. For example 
in [53] the FPL kernel modes has been derived using the same technique of Section 
3 in the case of the FPL equation. This originates the expression 

[3 L (l,m)= f dq*{\q\){[l^ 2 -[m^f}e^ m , 

JB(Q.2\tt) 

where I 1 ' and are defined as before. 

The equivalence of these representations can be proved by an indirect argument 
using the recent results in [29 . 
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iii) Using the symmetry properties of the Boltzmann kernel modes (jl5[) which are con- 
served in the grazing limit procedure, we obtain an almost complete set of sym- 
metries for the FPL kernel modes and that also the FPL kernel modes are real 
functions. 

4.3. Fast algorithms. We will show here that thanks to the particular structure of equa- 
tions (|2"4"|) and (|2HI22p . that define the FPL kernel modes and an approximation of the 
Boltzmann kernel modes in the grazing limit, the resulting spectral method can be evalu- 
ated with 0(N 3 log TV) operations instead of 0(N 6 ). 

We describe the method only for the FPL kernel modes, the extension to expressions of 
the form (I21j|22p will then follow easily. 

First we rewrite the resulting spectral scheme in the form, k — —N,...,N 

N 

9tfk= ^ f k - m f m @ L ( k - m i m )i 

m=-N 

where we assume that the Fourier coefficients are extended to zero for \kj\ > N, j — 1, 2, 3. 
Next we observe that the term m) splits as 

&L(!,m) = -V% / i qj «f(qy qm dq-k 2 [ \q\ 2 ^(q)e iqm dq 



(27) +7]%^/ ^(q) qj q h e^ m dq 

j,h=l JB(0,2\w) 

3 3 

:= kjFj(m) + k 2 G(m) + kjkhljh(m)> 
3=1 j,h=l 
where Ijh is a symmetric matrix. 

Thus the resulting scheme requires only the evaluation of convolution sums since it can 
be written as 

3 N N 

<kfk=/l kj fk-m f m Fj(m) + k 2 2J fk- m f m G(m) 

j m——N m=—N 
3 N 
+ ^ k 3 k h ^ fk-rn f m Ijh(m). 

j,h=l m.=-N 

It is well-known that transform methods enable to evaluate the previous expression in only 
0(N 3 log AT) operations [7], [16]. The computation of the terms Fj(m), G{m) and Ijh{m) 
can be done following the same lines presented in [24] . 

5. Uniform spectral accuracy of the method 

We explain in this section why the method in this paper is spectrally accurate, uniformly 
in the approximation of grazing collisions. 

Let us first introduce some notations. For any t > 0, /j\r(w, t) is a trigonometric polynomial 
of degree iV in v, i.e. f N (t) G P N where 

1P N = span {e lk ' v | - N < kj < N, j = 1, 2, 3} . 
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Moreover, let Vn ■ L 2 ([— n, 7r] 3 ) — > P N be the orthogonal projection upon P N in the inner 
product of i 2 ([-7r,vr] 3 ) (see (P]l) 

<f-P N f,<P>=0, V^elP N . 

We denote thei 2 ([-7r,7r] 3 )-norm by II/H2 = (< /, / >) 1/2 . In addition we define H£([-tt, tt] 3 ) 
where r > is an integer to be the subspace of the Sobolev space H t ([—tt,tt} 3 ), which con- 
sists of periodic functions [7] . 

Using these notations the method defined by equation (|13p can be written in equivalent 
form as 

(28) d t f N = Q x N (f N J N ) 

with the initial condition fw(v,t = 0) = Jo,n{v), where we have used Q x (f,f) instead of 
Qg(f, f) to denote the Boltzmann collision operator with relative velocity bounded by 2\ir. 
We point out that because of the periodicity assumption on /, and hence on Q x (f, /), the 
collision operator Q x {f, /) preserves in time the mass contained in the period. This can be 
proved directly using the property that f3(—m, m) = 0. On the contrary, momentum and 
energy are not preserved in time. 

Remark 5. Note that, the Boltzmann equation itself is non conservative in a bounded 
domain. Deterministic conservative schemes are usually constructed by modifying the col- 
lision mechanics or the numerical method in order to force the conservation properties in a 
bounded domain. However, it is clear that even with this modification the accuracy of the 
solution with respect to the original equation is guaranteed only if the size of the velocity 
domain is large enough. 

Following the same strategy as in [33], the spectral accuracy of the method will be a 
consequence of a uniform estimate on the true Boltzmann operator. This is the content of 
the following proposition. 

Proposition 6. Let B(q,0) = 3>(\q\)b(co$9), and 

Qb(/,<?) = / d«. / daB(v-v*,6)(g'J' '- /./). 
Js.3 Js 2 

Then 

(29) \\Q B (f,g)\\v> <C\\g\\ L 4f\\ H 2, 
where C is, up to a universal constant, 

r iv/2 



su P (i + M 2 )$(M) 



d9 sin6»6(cos6»)(l - cos0). 







Remark 6. 

i) By symmetrizing the kernel, we restricted the domain of 9 to [0, tt/2]. This is always 
possible thanks to the indiscernability of particles. 

ii) In the applications, we shall use this proposition with the support of / and g trun- 
cated over [—tt, 7r] 3 , so that the norm of g in L 1 is readily estimated by the norm in 
L 2 , and the constant C is finite even for $(|q|) = \q\ a , a > 0. 
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iii) In the singular case where = \q\ a , —3 < 7 < 0, similar estimates hold, if one 
allows to replace \\g\lL 1 by (say) \\g * This is estimated by a norm of g in 
some L p space, which in turn can be estimated (at least locally) by the norm of g 
in some Sobolev space H l . 

iv) This estimate is clearly, in the case which is of interest for us, uniform during the 
process of asymptotics of grazing collisions. 



Proof 

We use a duality argument : 



\\QB(f,g)\\ L i 



= SU P {/ 



QB{f,g)y, IMU2 < 1 



But 



J QB(.f,g)f = J ^dvdv* daB(v - v*,6)g*f(tp' - ip) 



< \\g\\ 



su P (i + M 2 )$(M) 



sup 

ii«6l 3 JR 3 



dv 



[ dab(cos6)f(<p' - ip) 
Js 2 



and 



/ dv iJ {V) 12 / da 5 ( cos °) & - <?) 
J R 3 l + \v-v*\ 2 J 



< \\f\\m 



1 



1 + \v - v*\ 2 



s 2 



dab(cos9)(ip' — ip) 



So it suffices to show that 
1 



J *\ JS 



da 6(cos6*)(<// — p) 



<cy\\ 



L 2 ■ 



H- 2 



By translational invariance, it suffices to prove this estimate for w» = 0. Since 1/(1 + \v\ 
is a smooth function of v, it suffices to prove that 



(30) 



<C\\<p\\h*, 



L 2 



where ip' = ip((v + \v\a)/2). 

We now note that (up to a multiplicative constant) 

f- pTt/2 f>1-R 

\ da 6(cos 9) (ip' — ip) = / sin#6(cos#) / d(j>(ip' — ip) 
Js 2 Jo Jo 



= / d9 sin 6»6(cos e)Vip(v) ■ (v' - v) d4> 
Jo 

+ de sin06(cos0) / / (1 - t)D 2 ip(v + t(v' - v))(v' - v, v' - v) dt d<j>. 
Jo Jo Jo 



Now, by symmetry, 



d(f>Vtp(v) ■ (v' — v) = V(p(v) ■ v(cose — 1), 
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S' 2 



da b(cos9)(ip' — ip) 



< C 



dO 



(cos0)(l - cos 6) 



\V<f{v)\\v\ + sup / dt(l-t)\D 2 <p{v + t(v' -v))\\v\' 

rr v>0 JO 

Then, (|30|) follows by noting that for each t E [0, 1], a E S 2 such that a ■ v > 0, one has 
/ dv \D 2 <p(v + t(v' - v))\ 2 < C J dv \D 2 ip\ 2 dv. (the Jacobian of v ^ v + t(v' - v) is bounded 
below). 



Next we state the consistency in the L 2 -norm for the approximation of the collision 
operator Q x (fJ) with Q^(/jv,/jv), 

Theorem 7. Let f E H 2 ([-n, tt] 3 ), thenar > 



(31) \\Q\fJ)-QN(fNjN)\\2<C [\\f-f N \\ H2 + 



where C depends on \\f\\% 



\\QHfNjN)\\H; 
N r 



Proof 

First, we can split the error in two parts 

\\Q X (f,f)-Q X N(fNj N )\\2 < \\Q X (fJ)-Q\fNjN)\\2 
+ \\Q X (fNjN)-Q X N (fNjN)\\2- 

Now clearly Q x (Jn, /jv) G and hence Q x (Jn, In) is periodic and infinitely smooth 



Q\h ■ / a ) - - Q-\ ( Wtf)||a < ^\\Q X (fNjN)\\ H; , Vr > 0. 



together with all its derivatives thus [7] 
(32) ur^rt- *-\ 

By application of Proposition [6] and of the identity 

Q A (/,/)-Q A (.9,.9)-Q A (/ + .9,/-3), 

we have 

\\Q X UJ)-Q X (fNjN)\\2 = \\Q X (f + fN,.f-fN)\\2<C\\f + f N \\l\\f-fN\\H} 

<2C 1 ||/|| 2 ||/-/ jv ||h|. 

■ 

Finally the following corollary states the uniform (in the grazing collision parameter) 
spectral accuracy of the approximation of the collision operator 



Corollary 7.1. Let f E H^([-%, tt] 3 ), r > 2 then 



(33) 



\\Q X (f,f)-Q X N (fN,fN)\\2 < 



C 



m + \\Q x (fNjN)\\ H r 
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Proof 

It is enough to observe that 

||/-/iv||^<^||/||^. 

■ 

Remark 7. From the previous corollary it follows 

I < Q\fJ)^> - < QUfr-Mv > I < ^ZalMla (ll/lk; + \\Q\fN, f N )\\ H; ) , 
and hence, by taking tp = v,v 2 , the spectral accuracy of the moments. 

6. Conclusions 

In this paper we have studied a Fourier spectral method that allows a numerical passage 
from the Boltzmann equation to the FPL equation in the grazing collision limit. In particular 
we have proved the uniform boundedness of the kernel modes for any collision operator 
with a singular kernel, provided the angular cross-section for momentum transfer © is 
bounded. This permits to show that in the small grazing collision limit the method provides 
a consistent discretization of the limiting FPL equation. The uniform spectral accuracy of 
the method with respect to the grazing collision parameter has been also given. Moreover, 
we have derived an approximated formula (|21[) that gives an intermediate asymptotic of the 
Boltzmann kernel modes which can be evaluated with fast algorithms similarly to the FPL 
modes studied in |24) . Numerical results based on this approximated formula and comparison 
with exact solutions are actually under development and will be presented elsewhere. 
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